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We investigate planar quantum electrodynamics (QED) with two degenerate staggered fermions 
in an external magnetic field on the lattice. We argue that in external magnetic fields there is 
dynamical generation of mass for two-dimensional massless Dirac fermions in the weak-coupling 
region. We extrapolate our lattice results to the quantum Hall effect in graphene. 
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I. INTRODUCTION 

Quantum electrodynamics (QED) in 2+1 dimensions 
is interesting as a model for several condensed matter sys- 
tems. In fact, quantum electrodynamics with two mass- 
less Dirac fermions could be relevant to describe the low- 
energy excitations of a single sheet of carbon atoms ar- 
ranged in a honeycomb structure called "graphene" (HQ. 
When graphene is immersed in a transverse magnetic 
field, the presence of Landau levels at zero energy leads 
to half-integer quantum Hall effect. Moreover, for very 
strong magnetic fields there is experimental evidence 
for the dynamical generation of a gap, which signals 
the spontaneous breaking of the chiral symmetry. In 
fact, it has been suggested that a magnetic field is a 
strong catalyst of chiral symmetry breaking in spinorial 
QED even at the weakest attractive interaction be- 
tween fermions. 

The aim of the present paper is to investigate, by means 
of non-perturbative Monte Carlo simulations, planar 
quantum electrodynamics (QED) with two degenerate 
staggered fermions in an external magnetic field. To 
make contact with the physical planar systems, we choose 
to work in the weak-coupling region. A preliminary ac- 
count of the results discussed in the present paper has 
been published in Ref. [5{. 

The plan of the paper is as follows. In Sect. HH for 
completeness, we briefly discuss our method to introduce 



background fields on the lattice and compare with differ- 
ent approaches in the literature. Section Hill is devoted to 
the discussion of our lattice Euclidean action. In Sect. II VI 
we present the results of our numerical simulations for 
two different values of the gauge coupling in the weak- 
coupling region. In Sect. [V] we extrapolate our results 
to the physical relevant case of the quantum Hall effect 
in graphene. Finally, our conclusions are relegated in 
Sect. [Vl] 



II. BACKGROUND FIELDS ON THE LATTICE 

The study of lattice gauge theories with an external 
background field has been pioneered in Ref. jf| 0] for the 
U(l) Higgs model in an external electromagnetic field. In 
the continuum a background field can be introduced by 
writing: 



A^x) -> A^x) + A^(x) . 



(1) 



In the lattice approach one deals with link variables 
U^{x). Accordingly, on the lattice Eq. ([T]) becomes: 



U^x) U,(x)U^(x) 



(2) 



where U° xt (x) is the lattice version of the background 
field A cxt (x). As a consequence the gauge action gets 
modified as: 
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S G [U] -»• S G [U] +5S[U,U C 



(3) 



where 6 S[U, U ext ] takes into account the influence of the 
external field 8-151. An alternative method, which is 
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equivalent in the continuum limit, is based on the ob- 
servation that an external background field can be intro- 
duced via an external current [16l - l20j : 

JT = d » ■ (4) 

The gauge action gets modified in an obvious manner: 



S 



B 



(5) 



where: 

S B = Jdx J™\x) A„{x) 

= -\ fdxFff(x)F vlx (x). (6) 

The background action Sb can be now easily discretized 
on the lattice. 

The main disadvantage of this approach resides on the 
fact that it cannot be extended to the case of non- Abelian 
gauge group in a gauge-invariant way. To overcome 
this problem, the background field on the lattice can 
be implemented by means of the gauge invariant lattice 
Schrodinger functional [2l|, [22| : 



VU e~ SG 



where the functional integration is extended over links 
on a lattice with the hypertorus geometry and satisfying 
the constraints (xt is the temporal coordinate) 



U k (x)\ xt=0 = U^(x) 



(8) 



We also impose that links at the spatial boundaries are 
fixed according to Eq. (JHJ. In the continuum this last 
condition amounts to the requirement that fluctuations 
over the background field vanish at infinity. 
The effects of dynamical fcrmions can be accounted for 
quite easily. In fact, when including dynamical fermions, 
the lattice Schrodinger functional in presence of a static 
external background gauge field becomes [23[ 



VUVtl)Vipe- {SG+SF) 



U k (L t ,x) = U k {0,x)=U^{x) 



VUe~ s ° det M , (9) 

IU k (L t ,S) = U k (0,3)=U^HS) 

where Sf is the fermionic action and M is the fermionic 
matrix. Notice that the fermionic fields are not con- 
strained and the integration constraint is only relative 
to the gauge fields. This leads to the appearance of the 
gauge invariant fermionic determinant after integration 
on the fermionic fields. As usual we impose on fermionic 
fields periodic boundary conditions in the spatial direc- 
tions and antiperiodic boundary conditions in the tem- 
poral direction. 



III. LATTICE PLANAR QED IN EXTERNAL 
MAGNETIC FIELD 



We are interested in planar quantum electrodynamics 
with Nf = 2 degenerate Dirac fields in an external con- 
stant magnetic field. As it is well known, Dirac fields are 



described non-perturbatively by the lattice Euclidean ac- 
tion using N flavours of staggered fermion fields \ 0] : 



JY 



S = S G + ^2^2xi(n)M n , m Xi(m) , 



(10) 



i— 1 n,m 



where Sq is the gauge field action and the fermion matrix 
is given by: 



M„ 



[U] = 
r) v (n) 



^ 2 

v=l,2,3 



{U v (n)8 mt n+i> - Ul{m)8 m ,n-i>} + TOo#m,n 



Vuin) = (-1)"- 



(11) 



where tuq is the bare fermion mass. Here we adopt the 
compact formulation for the electromagnetic field (for a 
detailed account see Ref. (25|). The gauge action is: 



S G [U]=P 



i 



l 



(fV(n) + [/t») 



(12) 



(7) where U^ v {n) is the plaquette and j3 



2 . The action 
Eq. (|10p with N = 1 flavours of staggered fermions corre- 
sponds to Nf = 2 flavours of 4-component Dirac fermions 

* m. 

To introduce an external magnetic field, we shall follow 
the lattice Schrodinger functional described in Sect. IIIII 
(for a different approach see Ref. [13]). Accordingly, in 
the functional integration over the lattice links we con- 
strain the spatial links belonging to the time slice Xt = 
to 



C4(£,a; t = O) = [/ fe cxt (20, k = 1,2 , 



(13) 



U£ xt (x) being the lattice version of the external contin- 
uum gauge potential. Since our background field does not 
vanish at infinity, we must also impose that, for each time 
slice Xt 7^ 0, spatial links exiting from sites belonging to 
the spatial boundaries are fixed according to Eq. (fT5|) . 
The continuum gauge potential giving rise to a constant 
magnetic field is given by: 



\x) 



3fe,2 



XiH , 



(14) 



so that: 
U? t (x) = l,U^ t (x) 



cos(eHxi) + i sm(eHxi) . (15) 



Since our lattice has the topology of a torus, the magnetic 
field turns out to be quantized: 



rr 27T 

e-H = —n cxt , 



integer , 



(16) 



where L is the lattice size. We recall once more that the 
fermion fields are unconstrained and satisfy antiperiodic 
boundary conditions in the timelike direction and peri- 
odic boundary conditions in the spatial directions. 
Our numerical results were obtained by simulating the 
action Eq. (fit)]) on L 3 lattice using standard hybrid Monte 
Carlo algorithm. 
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IV. CHIRAL SYMMETRY BREAKING 

We are looking for the dynamical generation of a gap 
for massless fermions. This corresponds to a non-zero 
chiral condensate in the chiral limit. Our strat- 

egy is to measure the fermion condensate with a small 
bare fermion mass tuq and then perform the massless 
limit mo — > in presence of a constant external magnetic 
field. Our simulations have been performed in the weak- 
coupling region with two different values of the gauge 
coupling j3 — 2.0 and /3 = 2.5. In fact, in the weak- 
coupling region we expect that the effects of the Coulomb 
interactions could be neglected allowing to extrapolate 
our numerical results to physical planar systems. 
We have performed simulations on lattices with L = 16, 
24 and 0.005 m ^5 0.03 with different strengths of 
the external magnetic field labelled by the integer n cxt 
according to Eq. (|T6|) . For each parameter set, to al- 
low thermalization we discard 10000 sweeps for L = 16 
and 7000 sweeps for L = 24. We collect about 50000 
hybrid Monte Carlo trajectories. To optimize the per- 
formance of the hybrid Monte Carlo algorithm, we tuned 
the simulation parameters to give an acceptance of about 
80%. The chiral condensate ('I"!') was estimated by the 
stochastic source method. In order to reduce autocorre- 
lation effects, measurements were taken every 10 steps 
for L = 16 and every 5 steps for L — 24 . Data were an- 
alyzed by the jackknife method combined with binning. 
In Fig. [T] we display the chiral condensate for different 
values of the lattice size, bare fermion mass, and mag- 
netic field strength for (3 = 2.0 . Note that, according 
to Eq. (fT6"|) . the strength of the external magnetic field 
depends on n cxt as well on the lattice size L. To avoid 
lattice discretization and finite volume effects, we have 
fixed the magnetic field strength such that the magnetic 
length satisfies the bounds: 



The best fit of the data to Eq. ([18)) in the scaling region 
gives: 




(17) 

We expect that in the continuum limit the relevant scale 
is set by the magnetic length. This means that the 

rescaled chiral condensate ^fjp would depend only on 

the scaling variable x = T "° . Actually, from Fig. [TJ 



where we display the rescaled chiral condensate versus 
the dimensionless scaling variable x, we see that in the 
region x > 0.05 data are rather scattered. However, in 
the region x < 0.05 our data seem to collapse to an uni- 
versal curve. This means that in this region, that we shall 
call the scaling region, the rescaled chiral condensate de- 
pends only on the scaling variable x. This allows us to 
extract the chiral condensate in the chiral limit mo — > 0, 
which corresponds to x — > 0, for a fixed strength on the 
external magnetic field. In fact, we try to fit the data in 
the scaling region < x < 0.045 according to: 



eH 
2tt 



CIq + di X 



m 



«0 

Xd.o.f. 



0.04399 ± 0.00131 
747 . 



ai = 9.759 ± 0.055 



(19) 

We note, however, that there are sizable violations of 
our scaling law as implied by the huge reduced chi- 
square. We believe that these scaling violations are 
mainly due to the fermion interactions with the electro- 
magnetic field, which could introduce a spurious depen- 
dence of the scaled chiral condensate on the dimension- 

2 

less ratio % . To check this point, we have performed 



numerical simulations by increasing the gauge coupling /3 
(which corresponds to a smaller e 2 ). In fact, in Fig. [5] we 
display the results of our simulations for f3 = 2.5. Again 
we see that the data for the rescaled chiral condensate 
seem to collapse to to an universal curve in the scaling 
region x < 0.05. Moreover, comparing Fig. [5] with Fig.[T] 
it is evident that the scaling violation are greatly reduced 
allowing a better extrapolation to the chiral limit. Fit- 
ting the data to Eq. ([18]) we find: 



a = 0.03544 ± 0.00084 



oi = 8.382 ± 0.032 



Xd.o.f. 



466 . 



(18) 



(20) 

Even though the reduced chi-square is quite large, we 
believe that our results are robust enough to allow the 
extrapolation of the chiral condensate to the chiral limit. 
As a consequence, we conclude that in the chiral limit 
the external magnetic field does induce a non-zero chiral 
condensate. From Eqs. (fT5)) and (|2"0"|) we find for the 
chiral condensate in the massless limit: 

— eH 

= — (0.03544 ± 0.00084) . (21) 

27T 

The non-zero value of the chiral condensate can be in- 
terpreted as the generation of a dynamical fermion mass 
which, in principle, can be extracted from the non-zero 
chiral condensate in the chiral limit. 

In the determination of the value of the chiral conden- 
sate, as given in Eq. (|21|) , we neglected a possible contri- 
bution present even in absence of the magnetic field. In- 
deed, in Ref. [28| it was shown that the chiral condensate 
is non-zero in the weak-coupling regime of compact pla- 
nar QED even at zero external magnetic field. In order to 
check the possible impact of this zero-field contribution 
on our determination of the chiral condensate, we observe 
that in Ref. [H[ two of us found /3 2 (f*) » 1.5X RT 3 for 
H = on a lattice with L = 12. This result implies, for 
(3 = 2.5, that = 0.00024, in lattice units. In this 

work the smallest value of the chiral condensate induced 
by an external magnetic field is obtained for n cxt = 1 and 
L = 16, which implies eH/(2ir) — 0.0625 and therefore, 
through Eq. (SI]), = 0.002215. The latter value is 

one order of magnitude bigger than the former and can- 
not be attributed to finite size effects, in consideration 
of the similar lattice sizes adopted in the two determi- 
nations. This allows us to safely neglect the zero-field 
contribution to the chiral condensate. 
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FIG. 1. Scaled chiral condensate versus the scaling variable x = 
the data Eqs. (fT8|) and (fl9]> in the scaling region < x < 0.045. 



y/eH/2n 



for fl = 2.0. The continuum line is the linear fit of 



V. EXTRAPOLATION TO GRAPHENE 

In this Section we attempt to apply our numerical de- 
termination of the chiral condensate in the chiral limit to 
graphene immersed in a transverse magnetic field. For 
the reader's convenience, we briefly discuss the remark- 
able quantum Hall effect in graphene. 
As is well known, graphene is a flat monolayer of carbon 
atoms tightly packed in a two dimensional honeycomb 
lattice consisting of two interpenetrating triangular sub- 
lattices (for a review, see Ref. [29]). Indeed, the structure 
of graphene has attracted considerable attention since the 
low-energy excitations are given by two Pauli spinors W± 
which satisfy the massless two-dimensional Dirac equa- 
tion with the speed of light replaced by the Fermi velocity 
vp ~ 1.010 s cm/s. The Pauli spinors can be combined 



into a single Dirac spinor VP = 




. Taking into ac- 



count the real spin degeneracy, we see that the low-energy 
dynamics of graphene can be accounted for by Nf = 2 
massless Dirac fields [10, HH . 

When graphene is immersed in a transverse magnetic 
field, the relativistic massless dispersion of the electronic 
wave functions results in non-equidistant Landau levels 



e n = sign{n)^j2\n\h^-eH , rc = , ± 1 , ± 2 ... (22) 

where eH > 0, e being the elementary charge (see 
Fig. |31 left). The presence of anomalous Landau lev- 
els at zero energy, £o — 0, leads to half-integer quan- 
tum Hall effect corresponding to quantized filling factor 
v = ±2, ±6, ±10 .... 

Recent studies of quantum Hall effect in graphene in very 
strong magnetic field H > 20 T (1 T = 10 4 gauss) have 
revealed new quantum Hall states corresponding to fill- 
ing factor v = , ±1 , ±4 [H, H3|- The new plateaus at 
v = , ±4 can be explained by Zeeman spin splitting. 
On the other hand the v = ±1 plateaus are associated 
with the spontaneous breaking of the symmetry in the 
n = Landau levels (the so-called valley symmetry). 
Indeed, these states are naturally explained if there is 
dynamical generation of a gap A (see Fig. [31 right). 
The gap Ao can be extracted from the measured activa- 
tion energy. In fact, in Fig. [4] we display the measured 
activation energy gap AE(i/ = 1) as a function of the 
magnetic field for the v = 1 quantum Hall states [34j . 
To extract Aq from the activation energy data, we need 
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the data Eqs. (Ti"8|) and (fl9]) in the scaling region < x < 0.045. 



y/eH/2n 



for fl = 2.5. The continuum line is the linear fit of 



to take care of the Zeeman energy which for strong mag- 
netic fields is no more negligible. To this end, we may fit 
the data to: 



AE(v = 1) = 2 



(MB) 



fJ 



H 



(23) 



where \xb is the Bohr magneton, g = 2 and Ao(-ff) ~ 
VH [3J]. Figure S] shows that, indeed, our Eq. (f23|) gives 
an excellent fit to the data. We find: 



M H ) = ( 13 - 57 ± 0.28) Kfc B y/H(T) , 



(24) 



where H{T) means that the magnetic field is measured 
in Tesla. 

Our strategy is, now, to relate the gap Ao to the chi- 
ral condensate. After that, using our determination of 
the chiral condensate on the lattice, we will estimate the 
gap and compare with the experimental determination 
Eq. 424]). 

To this purpose we follow Ref. (35|, where the hypoth- 
esis of rearrangement of the Dirac sea of graphene in 
an external magnetic field was used and the electron- 
electron Coulomb interactions were neglected. Note that 
in graphene the electron-electron Coulomb interaction, 



e 2 /r, in general is not small, so that this approximation 
could be questionable. A direct calculation gives (35j |: 

,- , . HceH 1 

= -2 A c 



2tt 



\/7T I 



A 



(25) 



where T(x) is the Euler gamma function and y) is 
the generalized Riemann Zeta function. For small gap, 
we may expand to the first order in Ao- Using r(i) = \pK 
and £(x, 1) = we get: 



HceH 
2-k 



2An 



2H^eH 



C(i). 



(26) 



This last equation relates the gap A to the rescaled di- 
mensionless chiral condensate. Using our determination 
on the lattice for the rescaled chiral condensate, we ob- 
tain: 



C( 



2vr 



ao 



(27) 
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FIG. 3. Schematic spectrum of Landau levels of graphene in applied magnetic field (left). Landau levels with dynamical 
generation of a gap Ao (right). The Fermi level is at sf = 0. 



where ao is given in Eq. (|20p . Finally, with the experi- 
mental value for the Fermi velocity we get: 

Ao(-ff) ~ 2.6 Kk B y/H(T). (28) 

Comparing Eq. ([28]) with Eq. ([24]) . we see that our esti- 
mate of the gap is about a factor five smaller than the 
experimental data. However, it is remarkable that we are 
able to reproduce the dependence on the external mag- 
netic field A ~ y/~H. 

VI. CONCLUSIONS 

We investigated planar quantum electrodynamics with 
two degenerate staggered fermions in an external mag- 
netic field on the lattice. Our numerical results seem to 
indicate that in an external magnetic field there is a non- 
zero chiral condensate in the chiral limit pointing to a 
dynamical generation of mass for two-dimensional mass- 
less Dirac fermions. 

We performed our simulations in the weak-coupling 
regime of the compact formulation of the lattice gauge 
action. As discussed in Ref. [28| . the non- compact for- 
mulation of the theory could have a different continuum 



limit than the compact one, the signature of this be- 
ing the different magnetic monopole dynamics, which in 
compact QED leads to an enhanced chiral condensate. 
As a matter of fact, in the compact theory the chiral 
condensate is non-zero in the strong-coupling regime and 
undergoes a crossover to a non-zero value in the weak- 
coupling regime, while in the weak-coupling regime of the 
non-compact theory it is compatible with zero. Although 
we believe that the numerical impact of this possible dif- 
ferent behaviour in the continuum should be negligible 
to our purposes, we plan to explicitly check this point by 
performing numerical simulations with the non-compact 
lattice gauge action. This will also permit us to make a 
comparison with the results of Ref. [2?} , where a differ- 
ent approach was adopted to introduce the background 
magnetic field on the lattice. 

We also tried to extrapolate our lattice results to the 
quantum Hall effect in graphene, since the low energy dy- 
namics of graphene is described by Nf = 2 massless Dirac 
fermions. Our non-perturbative Monte Carlo simulations 
allowed to confirm the dynamical breaking of the valley 
symmetry in the lowest Landau levels. Moreover, even 
though we greatly underestimate the dynamical gap, we 
were able to reproduce the dependence of the dynamical 
gap on the strength of the external magnetic field. 
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